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lon channels are transmembrane proteins that regulate and maintain the ionic concentrations across 
the cell membrane. Modeling the atomistic-level ionic flux through these channels is crucial for 
the understandingof several neurological diseases and related pharmaceutical discoveries. Exper- 
imental techniques now provide information about the channel’s physical structure which helps in 
developing realisticion transport models. lons entering a channel follow different trajectories as they 
traverse the channel; each associated with a certain probability. Quantities that explain these trajec- 
tories are the translocation and return probabilities, average lifetime, and spectral density (an exper- 
imentally accessible parameter) of ion number fluctuations. Theoretical analysis of ion transport has 
been limited to low-resolution continuum diffusion-based or kinetic-based models. Such analytical 
models fail to include key factors affecting the ionic conduction. In this paper, we extend previ- 
ous models by an electro-diffusion model incorporating the effects of electric field, energy barrier, 
and rate-limited association/dissociation of ions with surface charges inside the channel. Survival 


probability and spectral density are derived from the analytical model. 
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1. INTRODUCTION 


Ion channels are pores in a cell membrane that routinely 
allow the exchange of specific ions across them. These 
channels play a central role in regulating the chemical 
and electrostatic gradients across the cell membrane which 
guides the movement of ions in/out of the cell. The physics 
of ion movement through the ion channels (when open) is 
simple: electrodiffusion and convection. The challenge lies 
in predicting the transport of ions of different types and 
concentrations across ion channel structures that can have 
unique topographies. Additionally, channel proteins have a 
fixed charge that demands counter ions in a nearby ionic 
atmosphere. Imbalances in ion concentrations create strong 
electric fields, and the non-equilibrium behavior of ions in 
such fields is crucial for regulating the open/closed states 
of the channel. Recent experimental advances and new 
molecular-level simulation modelsare helping us under- 
stand the complex mechanisms involved in ion transport 
through channels.! 

Till recently, theoretical analysis of ion transport has 
been limited to diffusion-induced transport of metabolites.” 
Even though diffusion remains to be the main transport 
mechanism, electric field-induced drift velocity adds to 
the momentum of the transporting ions. This transmem- 
brane potential gradient, which drives ion channel gating 
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and permeation, should be considered an integral player 
in any ion channel modeling analysis. Besides drift and 
diffusion, ions are also being randomly associated or dis- 
sociated with the charged surface proteins in the channel’s 
vestibule. These surface charges are known to influence 
the gating, conductance, and toxin-binding effects of KcsA 
ion-channels, and their actual physical locations are pre- 
sumed to be in the selectivity filter.> So, it is also impor- 
tant to include the effects of association/dissociation of the 
metabolites in any analytical model which affects the mea- 
sured spectral density of ion number fluctuations in an ion 
channel. 

An ion entering a channel can either return to the source 
reservoir or traverse the channel’s vestibule and reach the 
sink reservoir successfully. Key quantities that describe the 
process would be the translocation and return probabilities, 
average translocation lifetime inside the channel and the 
spectral density of ion number fluctuations during the ion 
flow.‘ In this work, a physics-based model is presented that 
incorporates the effects of ion mobility, constant electric 
field, and association/dissociation of ions inside a cylindri- 
cal ion-channel. The main focus is on ion transport through 
a channel, including the interaction of an ion with sur- 
face protein charges but ignoring ion—ion interactions. In 
the next section, we derive a general expression of the 
probability density function (p.d.f.) related to the particle 
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number fluctuations inside an ion-channel. The p.d-f. is 
derived in the axial direction and incorporates rate con- 
stants related to the association/dissociation of ions. The 
generalized expression for the p.d.f. is represented in terms 
of Green function. In the following section, we derive 
the exact Green functions for diffusion-induced ion trans- 
port in the presence of a constant electric field. Using a 
transformation, we demonstrate a method of handling the 
effects of carrier mobility and electric field in such mod- 
els based on Green functions. Lastly, we derive the exact 
expression of the p.d.f., survival probability, and spectral 
density of particle number fluctuations. The results are 
plotted using Maple™ and simulation plots are provided. 
Such models can be extended to describe other charge 
transport processes in nature (besides ion transport) such 
as charge migration along polymers or translocation of 
DNA chains through artificial nanopores. 


2. FORMULATION 


Consider acylindrical ion channel of length L with two 
bulk reservoirs: source and sink (Fig. 1). The ions flow 
from source to sink guided by a concentration gradient and 
an electric field. We assume that the ions are very few in 
number inside the channel and mutually non-interacting. 
The ions can, however, interact with a single surface site 
inside the channel through association/dissociation. Let us 
denote by p(x, 1), the probability density of a free ion 
at an axial distance x along the channel and at a time f, 
whereas q(x, t) denotes the probability density of finding 
an immobile ion bound temporarily with surface protein 
charges. These functions obey the following equations,°* 


Op(x, t) 
Ot 


= £p(x, t)—k,6(x — a) p(x, t) 
+k ,6(x — a)q(x, t) (1a) 


dq(x, t) 
ot 
where x = a is the contact site of association/dissociation, 


k, is the rate of association, and k, is the rate of dis- 
sociation. In terms of the potential energy U(x) and 


= k,6(x —a)p(x, t)—k,6(x—a)q(x, t) (1b) 


: bound ion 
free ion 
source @ ® sink 
reservoir @ reservoir 
ay surface 
7 charge 


xX 


Fig. 1. A sketch of the cylindrical ion channel of length L and two 
reservoirs at either end. An ion traverses from the source to the sink 
reservoir. The free ions can occasionally bind reversibly with surface 
protein charges embedded in the channel’s wall. 
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the diffusion constant Dp, diffusion operator £ in one- 
dimension is given as,’ 
0 0 

£= Dee se (2) 
To obtain solutions of the two functions, p(x,t) and 
q(x, t), we introduce the Green function G(x, s|x,) in the 
Laplace domain satisfying (s —£)G(x,5|x,) = 6(x—a), 
with an initial impulse Q,5(x—x,) of amplitude Q, at 
an axial distance x = x,. With this initial condition, we 
assume an ion impulse is externally applied precisely at 
x =x, and time progression of this impulse is analyzed 
to obtain the system impulse response. Furthermore, an 
absorbing boundary condition for both source and sink 
reservoirs is assumed. As shown in Appendix A, the gen- 
eralized solutions of p(x, ft) and q(x, ft) in the Laplace 
domain can be written as: 


2 G(a,s | x,) sk, 
PA $150) = Tek T(t )G( 512) an) 
x G(x, 5 | a)+G(x,5|x,) (3a) 
: G(a,s | x,) 
9% 51%) = Te +k Gs 1a) 
x (= )oe-0 (3b) 


where p(x, s | x,) and qg(x,s5|x,) denote the Laplace 
transforms of p(x, ft) and q(x, t) respectively with an ini- 
tial impulse at x = x,. Equations (3a) and (3b) describe 
the probability density functions (p.d.f.) corresponding 
to an ion in the unbound state or in the bound state. 
Next, we will derive the expression of the Green func- 
tion G(x, s|.x,) for substitution into Eq. (3a) to obtain 
the exact analytical expression forthe p.d-f. of a free ion, 
P(x, s|x,). This result willthen be used to derive the spec- 
tral density of ion number fluctuations in the ion channel. 


3. GREEN FUNCTIONS 


We assume, for simplicity, one-dimensional transport of 
ions in a constant electric field ¢. The current density and 
the continuity equations can be written as,? 


Jp = qbpeg (x, t) — qDpdg(x, t)/dx (4a) 
OJp/Ox + Gdg(x, t)/dt = —Gg(x, t)/Tp (4b) 


where Jp is the current density, g is the ionic charge, Dp 
is the diffusion constant, wp is the mobility, Tp is the 
recombination lifetime, ¢ is the constant electric field, and 
g(x, t) is the Green function in the time-domain. We are 
interested in finding an expression for G(x, s | x,) defined 
in the previous section, which is the Laplace transform of 
g(x, t). Substitution of Eqs. (4a) in Eq. (4b) gives, 


po 8b) dg(x,t) g(x,t) _ dg(x, t) 
PT ayz Bre ay «oF 


(5) 
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We now introduce the transformation, g(x,t) = 8E+6 
g' (x, t)exp[—t/7p], that helps us to rewrite Eq. (5) as: 7E46 4 
0g! (x,t Og’ (x,t Og'(x,t 6E+6 5 2 200s 
D, g(x 0) ge ) _ 9g'(x, 0) (6) ane 
Ox? Ox ot 5E+6 J c 0.07ns 
= d 0.02ns 
As shown in Appendix B, the two solutions of Eq. (6) are: - 4E+6 4 e 0.01ns 
3E+6 | 
G(x, s|x<x,) 
arse, 2E+6 | 
PoSinh,/'s/Dp(L—x)+ Qo,/s/Dpsinh,/'s/Dp(L—xo)sinh,/s/Dpx 
= s-sinh,/s/DpL IE #6) 


(7a) 
G(x,s|x>x,) 


Posinh,/s/Dp(L—x)+ Qoy/'s/Dpsinh/s/Dpxp sinh /'s/Dp(L— x) 
_ s-sinh /s/DpL 


(7b) 


where G(x, 5|x <.x,) and G(x,s5|x >-x,) together repre- 
sent G(x, 5|x,) which is the Laplace transform of g(x, f), 
Po is the bulk concentration of the source reservoir, and 
Qo is the amplitude of the initial pulse. It is worth men- 
tioning that in Eqs. (7a) and (7b), the effect of mobility 
and electric field are not directly obvious because of a 
(x’, t') space-time transformation made during the deriva- 
tion (see Appendix B). However, after Laplace inversion, 
we can revert to the original t(x, ft) space-time domain 
where the terms containing mobility and electric field are 
obvious. The inverse Laplace transform of Eq. (7) gives 
the Green function solutions in the time domain without 
(Eq. (8a)) and with (Eq. (8b)) the effect of the electric field 
as follows: 


eB Det! 


ye sin (Bx') 


n=1 n 


a: 
a(x’,t) = po} 1-= 
7 


2 ioe) —B°Dpt' 
205+! np ain (ie) (8a) 
n=1 n 

where 8 = n7r/L. Substituting back the initial transforma- 
tions from (x’, t’) to (x, f), 


2 2 Ee BDet 


rf (Bs— rer)| 


n 


20,60 sa 


TL 
x sin (BXo) (8a) 


In Figure 2, we plot the time variation of the Green func- 
tion solution g’(x’, t') from Eq. 8(a). The channel length 
is assumed to be 40 A and the initial pulse, 5(x — x9), 
is applied in the middle of the channel (xg = 20 A). As 
shown in the figure, the Green function spreads uniformly 


sin (B(x —mpet)) 
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axial distance (A) 


Fig. 2.. The figure shows the diffusion-influenced time variation of the 
Green function from Eq. (8a). The channel was assumed to be 40 A long 
and the initial pulse is applied in the middle of the channel. The electric 
field is not applied here and the pulse decays with time as a result of 
diffusion. 


to either side of xg, gradually diminishing in peak value 
due to diffusion. In Figure 3, we plot the same time varia- 
tion of the Green function solution g(x, t) from Eq. (8b). 
The channel is 40 A long and the initial pulse is applied 
in the middle of the channel. In this case, we notice that 
the Green function slowly migrates to the right end of the 
channel caused by the right-side directed electric field and 
spreads out due to diffusion. 

To show the effect of electric field and surface sites on 
ion transport, we simulated a cylindrical channel (40 A 
long and 6 A wide) using a SILVACO simulation soft- 
ware. In this solid-state simulation tool, silicon is used to 
mimic the conducting water continuum and silicon oxide 
to mimic the non-conducting protein walls. The simula- 
tion parameters (carrier concentration, ionization degree, 
and diffusion coefficients) are adjusted in each conducting 
region to emulate the realistic motion of potassium ions 


0 5 10 15 20 25 30 35 40 
axial distance (A) 


Fig. 3. The figure shows the time diffusion and electric field-influenced 
variation of the Green function from Eq. (8b). The channel was assumed 
to be 40 A long and the initial pulse is applied in the middle of the 
channel. The pulse decays with time as a result of diffusion and also 
propagates to the right due to the applied electric field. 
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inside the channel.!° The simulations are performed using 


the discretized PNP Eqs. (4a) and (4b). Figure 4 shows 
decreased ion concentration inside the cylindrical chan- 
nel, as predicted by Brownian dynamics.!° Figure 5 plots 
the potential profile in the presence of an electric field 
and increasing surface charges. As shown in the figure, 
an energy barrier inside the channel is transformed to a 
potential well by the inclusion of surface charges, thus 
increasing the chances of ion permeation. The electric field 
further lowers the potential well towards the sink, making 
it easier for ions to traverse the channel region. 


4. SPECTRAL DENSITY 


We now substitute the Green function G(x,5 |x < x,) 
and G(x, s|x > x,) from the previous section (Eqs. (7a) 
and (7b)) into Eqs. (3a) and (3b) to obtain the exact 
p.d.f. of a free ion, p(x, s|x,) and that of a bound ion, 
q(x, 5|x,). Assuming a perfect sink on either side of the 
channel (pp = 0) and a unit impulse as the initial state 
(Q, = 1) described in Section 2, the p.d.f. of a free ion 
P(x, 5| x,) are given by: 


Pi (x, 5| x <x,) 


= (Gro > sinh (/ ects) 
x sinh (/=)) /(s| <0 +4 sin (=) 
+ hay | sinh (5-449) sinh Je) |) 


(9a) 
P2(x,5| x > X,) 
S. S 
= | (s+k,),/ —sinh{ ./—(—-L+ x) 
Dp Dp 
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Fig. 4. Comparison of the potassium ion concentration in a cylindri- 
cal channel (40 A long and 6 A wide) obtained from this model and 
Brownian dynamics!” in the absence of electric fields. 
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Fig. 5. In the presence of an electric field and increasing surface charge 
concentrations, the energy barrier is lowered and it is easier for ions to 
pass through the channel region. 


com) [CL -wssoa (fi) 
#ky| > sinh ( RCE+m)) sit (55-0) ]) 


(9b) 
The correlation function of interest is the survival proba- 


bility which is given by, 


cilad=faicatlayt f pxletlx) — (10) 


The survival probability is an indicationof how successful 
an ion is to traverse the ion-channel from the source to the 
sink reservoir (Fig. 1). In other words, survival probability 
is the probability that a free ion will remain free (or dis- 
sociated) and not associate with surface charges on the 
channel wall. In Laplace domain, it can be expressed as: 


L{C(|x,)}=C(s|x,) =f Cr lxJe“ae (11) 
0 
Using Eqs. (9)-(11), the survival probability is obtained as: 


C,(s | x0) 


- ( +k) — sinh (5-4) + sinh (ec se *)) 
san (Ee) floss (ED 
affects) 


(12) 


In a typical electrophysiological experiment, one of the 
key measured quantities is the power spectral density, 
S(f), of the ion number fluctuations through the chan- 
nel. As reported through experimental reports, this noise 
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induced by ion transport exceeds the background noise by 
orders of magnitude.!° As such, the spectral density, S(f), 
is a valuable source of information about ion movement 
in a channel’s vestibule.'':!* Fortunately, we can relate the 
spectral density to the survival probability as: 


S(f) = 4Re{C(s = 27fi | x,)} (13) 


As shown above, spectral density is the frequency-domain 
Fourier representation of the correlation function or sur- 
vival probability in time-domain. For simplicity, let us 
assume that the initial impulse is given at the middle of the 
channel (x, = L/2). Using Eqs. (12) and (13), the expres- 
sion for spectral density is obtained as: 


Sif) = 4.Re| (202mfi+,) 


. [cosh (5v2afi70,)1) Z i}) 
i (enn 2008 +k,) cosh (5(v2mh7D,)L) 


+ (y 2fi7Dp)sinh (5 (/2mfi7Dp)4)})} (4) 


To demonstrate some applicability of the above equation 
and its dependence on the rate constants, we chose a set 
of rate constants and calculated the spectral density. The 
channel dimensions are kept same as before. Figure 6 plots 
the spectral density from Eq. (14) for the case of a fixed 
k, = 10-3 and values of varying kj. Figure 7 plots the 
spectral density of particle number fluctuations obtained 
from Eq. (14) with fixed k, = 10-> and varying k,. At 
low frequencies (f < 100 Hz), the spectral density is flat 
and then decays with increasing frequency. The frequency 
roll-off is typically represented as 1/f* where the exponent 
a can be between | and 2.” In both figures, the bottom 
curves correspond to the case of free diffusion where the 
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Fig. 6. The spectral density of particle number fluctuations obtained 
from Eq. (14) with fixed k, = 10* and values of k, varying from bottom 
to top: infinity (free diffusion), 9.91 x 10*, 9.97 x 10*, 9.93 x 10° and 
1.01 x 10°. 
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Fig. 7. The spectral density of particle number fluctuations obtained 
from Eq. (14) with fixed k, = 10? and values of k, varying from bottom 
to top: infinity (free diffusion), 9.935 x 10', 9.96 x 10° and 1.01*. 


rate constant of association or dissociation is assumed to 
be infinity. For finite values of the two rate constants, the 
low-frequency spectral density increases compared to the 
case of free diffusion. This suggests that free ions associate 
and dissociate with the binding sites with a finite rate, 
which increases the fluctuations in ion number and the 
spectral density. A comparison of both the figures shows 
that the roll-off slope is.more sensitive to the value of 
dissociation rate constant k, than that of the association 
rate constant k,. 

As mentioned earlier, the spectral density can be exper- 
imentally determined while other parameters such as sur- 
vival probability and p.d.f.s are difficult to determine 
directly from experiments. It is also not possible to deter- 
mine the rates of association (k,) and dissociation (k,) 
directly from experiments. However, it is possible to 
extract the above-mentioned parameters by fitting the pre- 
sented model with experimental data on current noise- 
generated spectral density of ion channels.'* The values of 
rate constants are chosen here to keep the noise magnitude 
(0.01 to 10 x 10-7’ A?/Hz) and frequency range (10 Hz to 
10 kHz) roughlyclose to those observed in actual ion chan- 
nel experiments.’? It is worth mentioning that the mag- 
nitude of this measured noise is at least 100-fold higher 
than shot noise or Johnson noise of the channel. Eq. (14) 
represents the general expression for plotting the spectral 
density which may be tailored to a different set of rate 
parameters or channel dimensions depending on the appro- 
priate ion channel under study. 


5. DISCUSSION 


Here, the modeling of the ion flux is inspired by the 
geminate excited-state reversible dissociation reaction that 
requires introduction of Green functions for depicting 
the transition probabilities of an ion’s association or 
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dissociation.2 Theso-called trapping problem, where a 
molecule diffuses among many static traps, has been stud- 
ied in the areas of chemical kinetics and material physics.® 
These studies, however, have been limited to solution- 
phase chemistry (where there are reversible diffusion- 
influenced reactions in the excited electronic state) or 
semiconductor physics (where various devices employ 
the electronic trapping/excitation for their working).!° 
In our models, we have extended the application of 
such diffusion-influenced reversible trapping problems to 
describe the ionic transport through ion channels and 
obtain the spectral density of ion number fluctuations. 

There are some limitationsinthe derived analytical 
model. In practical, diffusion coefficient used in our mod- 
els might not be a constant (in fact, they are time-varying). 
This is because the restricted diffusion of molecules in 
biological systems is affected by the presence of differ- 
ent types of barriers/obstacles. This tends to the slow- 
down of diffusion at long times as compared to barrier-free 
diffusion.'* In addition, our calculations assumea single 
ion channel is in the open state. Realistically, there may be 
multiple ion channels that may open and close for random 
time intervals. A single channel may undergo transitions 
between open, closed, and inactivated states. This chan- 
nel flickering may contribute to increased spectral density 
with several corner frequencies which may, however, be 
challenging to incorporate in the present model. 

With new analytical methods being. developed _ and 
increased computational capability, simulation models. of 
ion permeation are becoming increasing sophisticated 
andtypically fall into three categories: brownian dynam- 
ics, molecular dynamics, and continuum electrostatics.'!> '¢ 
Exact molecular dynamic solution of ion transport in a 
channel structure is ideally possible but is virtually impos- 
sible with today’s computational power. For example, to 
simulate | ps of real time for an ion channel model having 
1020 amino acid residues, 13,000 water molecules, and 
27 sodium atoms would require two hours of computation 
time. 

To reduce the computer time, continuum approxi- 
mation is selectively applied in these simulations to 
relatethe ion permeation characteristics to the physical 
channel structure.!® !7 With recent explosion in the amount 
of available genomic and structural data about mem- 
brane proteins, computationally efficient approaches are- 
being investigated to facilitate a greater understanding of 
channel properties and functioning.!? New computational 
approaches such as semi-Markov chain model!? and six- 
state kinetic model”° have been proposed to predict ion 
permeation with less computational complexity and com- 
parable accuracy to Brownian dynamics. Recently, syn- 
thetic nanopores have been used as experimental test-beds 
to validate several theories on ion channel transport such 
as preferential ion selectivity.”! In this respect, analytical 
models of ion transport could help complement available 
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simulation models and predict a number of key parame- 
ters. This work is an attempt to model the noise fluctua- 
tions during electrostatics-based ion transport in a channel 
environment and relate it to an experimentally measured 
parameter such as the spectral density. Future work is 
needed to incorporate the effects of ion selectivity, ion—ion 
interactions,hydration energy, and multichannel activity in 
such analytical models with regards to electrophysiological 
data.” 


6. CONCLUSIONS 


We presented an analytical method for calculating key 
parameters of ion transport through ion channels including 
survival probability of a free ion and spectral density of ion 
number fluctuations. Besides diffusion, the model includes 
the effects of electric field and reversible binding of free 
ions with the protein surface charges. The Green function 
is used to calculate the exact probability density functions 
of a free or bound ion within the channel’s vestibule. The 
general expression for spectral density shows the interplay 
of several key parameters governing ion transport, namely 
rates of association/dissociation, channel dimensions, ion 
mobility, and frequency. We saw that the spectral density 
has a 1I/f slope and increases with finite rates of associa- 
tion/dissociation between the ions and surface charges. 


Appendix A 


General Derivation of the p.d.f in terms of the Green 
Function 


Without a loss of generality, we can assume an initial state 
of free ions at a distance x,:° 


p(x, 0) = 6(x—x,), 
In Laplace domain, equations la and 1b can be written as: 
(s—£)p(x, 8 | x,) —8(x—x,) 
=[-k.P(4, 5 | xo) + kaq(a, 5 | x,)]6(x—a) (A2a) 
sq(x,5|X,) 
=[—kiq(a, 5 | x,) + kaP(a, 5 | x,)]6(x—a) (A2b) 


q(x,0) =0 (Al) 


where p(x, s | x,) and g(x, s|x,) denote the Laplace trans- 
forms of p(x, t) and q(x, ft), respectively with an initial 
condition given by Eq. (A1). We introduce the Green func- 
tion satisfying the equation: 


(s—£)G(x, 5 | x,) = 6(x—-a) (A3) 
Then, in terms of the Green function, Eq. (A2a) becomes: 


P(x, | Xo) = [—kap(a, 8 | x5) +kag(a, s | x,)] 
x G(x, 5|a)+G(x,5|7,) (A4) 
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From Eqs. (A2b) and (A4), we arrive at the following 
identities: 

G(a, s | x,) 
1+ (ska/(s+kq))G(a, 8 | a) 


P(a, 5 | x,) = q(a, 5 | x.) 


= aa P(a, 5 | x,) (A5a) 
—k,P(a,5|x,)+kaq(a, | x,) 
= G(a, 5 | x,) 
1+ (sk,/(s+k,))G(a, | a) 
Sk, 
x (-=:) (A5b) 


Using Eq. (A5b), the two probability density functions are 
given by Eqs. (3a) and (3b). 

Appendix B 

Derivation of the Green Function 

Rewriting Eq. (6): 


Ag(x.t) .08'(x.t) _ a8'(x, 1) 
Ax? ae at 


Dp (B1) 


We use a transformation x’ = x — ppet and t’= +t on 
Eq. (B1) yielding: 


09! (x', t') _ dg! (x, t') 


D 
Pax’? ar’ 


(B2) 
Using the Laplace transform on Eq. (B2) as: 
Lea's) = 808) = fa’. exp(—st')ar’ (B3) 


with the initial condition g’(x’, 0) = Q,6(x' — x9) gives: 


dg’ (x', s) Soi), Qo , 

ax? —-«D,® (x5) = Dee —Xo) (BA) 
The above initial impulse Q56(x—x,) has amplitude of 
Qo at an axial distance x = x, along the channel. Because 
of this generalized definition, one can vary the amplitude 
and location of the initial impulse according to specific 
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needs. The following are the boundary conditions: 


g(x’ =Xo, 5) lex, = g(x’ =Xo> 5) lox, 


a@(s's5)) ABS), _ Bo 
dx’ 70 dx’ 70 D> 


8 (x, Svar =0 8(%',5)lvo0=Po/s (BS) 


It was assumed that the source reservoir has an ionic- 
concentration po, the sink reservoir is perfect, and continu- 
ity is maintained before and after the site of initial impulse 
x =x,. We can now solve Eq. (B4) and obtain the Green 
function G(x, s5|x,) given by Eqs. (7a) and (7b). 
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